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EXPERIMENTS IN LOW ENERGY PLASMAS 


Microwave Studies of Electron Energy Loss Processes in Afterglow Plasmas 
Electromagnetic Wave Propagation through Magnetoplasmas 
Plasma Transport and Stability in Magnetic Fields 

E. H. Holt and J. H. Noon 
Rensselaer Polytechnic Institute, Troy, New York 


SUMMARY 


This report reviews research carried out in the 
Plasma Research Laboratory (renamed the Plasma Dynamics 
Laboratory in August, 1966) at Rensselaer Polytechnic 
Institute under NASA grant no. NsG ^t-8. The work is 
described under three major headings, microwave studies 
of electron energy loss processes in afterglow plasmas, 
electromagnetic wave propagation through magnetoplasmas 
and studies of plasma transport and stability in the 
presence of magnetic fields. 



INTRODUCTION 


The results are reported, of a comprehensive research program on the 
properties of collision dominated plasmas. Cross section data, transport 
coefficients (including the special case of time -dependent transport), 
propagation coefficients for electromagnetic waves and plasma oscillations 
are discussed. The current convective instability is studied in transverse 
fields and a strong stabilizing effect is reported. 

Understanding of the behavior of collision dominated plasmas is 
important in ionospheric work and in certain types of plasma accelerator. 
Because they are comparatively well understood such plasmas provide a 
favorable starting point for studies of plasma stability and containment 
in various magnetic field configurations. 

The work reported here includes a detailed study of the approach to 
thermodynamic equilibrium of a nitrogen afterglow plasma. Microwave inter- 
ferometry provides an account of the variation of the plasma density and 
microwave radiometry adds important information about the changing energy 
of the electron gas. Finally, a time resolved measurement of the velocity 
distribution of the electrons is made in order to provide a complete 
description of the electron behavior. 

A theoretical formulation is made of the effect of the wave polarization 
on the propagation of electromagnetic waves in magnetoplasmas. An experi- 
mental demonstration of the theory is provided using microwaves propagating 
in free space. The use of the measurement of wave polarizations as a 
diagnostic tool is enhanced by the use of guided waves. A plasma cell with 
turnstile junctions suitable for such studies is described. 

Collision dominated plasmas in strong magnetic fields are subject to 
the helical instability. Measurements of the oscillations which accompany 
this instability have resulted in a definitive understanding of the modes 
of the instability. A pre-excitation technique has provided a check of the 
linear theory which avoids the non-linear character of the full grown 
instability. The effects of transverse magnetic fields have been examined 
and it has been found that a minimum B field has a strong stabilizing 
effect. 

A microscopic theory of time -dependent transport phenomena is described. 
The advantages of this approach compared to modified macroscopic approach 
are discussed. 
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RADIATION TEMPERATURE 


MICROWAVE STUDIES OF ELECTRON ENERGY LOSS PROCESSES 
IN AFTERGLOW PLASMAS 


In order to monitor the electron temperature of the plasma as a function 
of time a gated microwave radiometer was developed (Ref. l). The radiometer 
has a time resolution of 1 microsecond and a temperature sensitivity of around 
10° Kelvin when measuring temperatures near room temperature. 

Two special waveguide cells were developed for microwave studies of after- 
glow plasmas (Refs. 2 and 3). These bakeable metal cells have advantages over 
glass discharge vessels and reduce uncertainties in the boundary conditions to 
be applied. 



TIME (/»••<:) 


Figure 1. Radiation temperatures 
measured for different gas pressures 
and different energy input conditions 
in the active discharge. Pulse widths 
from 3 microseconds to 20 microseconds 
and discharge currents from 0.5 to 
17- 5 amps were used. 
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Radiation Temperature Measurements 


Using the gated microwave radiometer, noise radiation from the electrons 
of the plasma was monitored as a function of time both in the active discharge 
and in the afterglow period following termination of the discharge. Radiation 
temperatures were obtained by comparison with a standard noise source, as 
described in Reference 1. 

For the nitrogen afterglow the initial fall. in the radiation temperature 
after termination of the active discharge agrees with what would be expected due 
to electron collisions with the neutral molecules. However, the radiation temper- 
ature then rises with time, reaching a maximum and then falling off slowly. 

Typical results for different pressures and different discharge conditions are 
shown in Figure 1. It is obvious that a source of energetic electrons must exist 
in the afterglow and that the energy input in the active discharge determines the 
magnitude of the effect. This work has been reported as a RASA Technical Note 
(Ref. k) and as an article in the Physical Review (Ref. 5). 

The effect of a small proportion of CC >2 on the radiation temperature of 
the nitrogen plasma is quite striking, see Figure 2. There is a near resonant 
energy match between the v = 1 vibrational level of the Ng molecule and the 
(00° l) rotational level of the COg molecule leading to strong quenching of the 
N*o ( v = l) level (Refs. 9 and 10). Thus, it is likely that excited molecules 
play an important role in the afterglow. 



Figure 2. Changes in radiation 
temperature of the nitrogen 
afterglow caused by admixture 

of a small percentage of C0 2 
in the discharge. 
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Electron Collision Probability- 


Calculations by other workers (Kefs. 11 and 12) of the probability for 
momentum transfer P c of low energy electrons ( ^ leV) in nitrogen as a 
function of energy, based on microwave conductivity measurements, have not 
included direct measurement of the electron temperature in the plasma. Results 
obtained by waveguide transmission (Ref. ll) and resonant cavity (Ref. 12) 
techniques differ by a factor of five at low energies. Although the measure- 
ments were carried out at different times in the afterglow, both sets of cal- 
culations included the assumption of a Maxwellian distribution of electron 
velocities whether or not the mean electron energy is greater than that of 
the neutrals. 

Our data indicates that calculations of P Q based on the assumption 
that only 5° microseconds are required for the electrons in the nitrogen after- 
glow plasma to reach the temperature of the neutral gas is questionable. This 
point was raised by Formato and Gilardini (Ref. 13) who observed a slow relaxation 
of electron temperature in nitrogen but did not carry out measurements at times 
close to the termination of the active discharge. 



Figure 3- Ratio of real and 
imaginary parts of the micro- 
wave conductivity ( a «• ) 

plotted as a function of radi- 
ation temperature for post-dis- 
charge times ranging from 30 to 
65 microseconds. This should 
fall off monotonically but does 
not do so, indicating that the 
calculations based on an assumed 
Maxwellian distribution are 
incorrect. 
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We have calculated P c as a function of energy "by equating the measured 
radiation temperature to the electron temperature . Although the results were 
generally close to those derived from transport coefficients (Ref. l4) and 
resonant cavity measurements (Ref. 12) they showed a wide scatter. Furthermore 
investigation of the behavior of the conductivity ratio <r / er (Figure 3), 
reveals that it is a multi-valued function of the radiation temperature instead 
of decreasing monotonically. A calculation shows that electron- ion collisions 
cannot explain the effect at these temperatures (approximately 2000°K). The 
most likely explanation is that the radiation temperature cannot be equated 
to the electron temperature and that the form of the electron velocity distri- 
bution is not Maxwellian in the early post-discharge period, as assumed in the 
calculations. This work has been reported in a NASA Technical Note (Ref. 15 ) 
and a Physical Review article (Ref. 5)- 


Non- Maxwellian Electron Velocity Distribution 


We have applied a microwave technique first described by Fields et. al. 

(Ref. l6) to determine the electron velocity distribution in the nitrogen after- 
glow plasma from the noise emission spectrum. If a longitudinal magnetic field 
is applied, measurement of the noise power from the plasma, as a function of 
magnetic field, for a range of field strength such that the electron cyclotron 
frequency is close to the microwave observational frequency, may be analyzed to 
determine departures from a Maxwellian distribution. A magnetic field of around 
3,250 gauss is required for measurements at X-band frequencies. The experimental 
results obtained reveal that a non-Maxwellian^situation exists in the nitrogen 
afterglow for periods of time of more than lCr electron- neutral mean free times 
and that the form of the distribution function changes rapidly following termi- 
nation of the active discharge. 


The radiation temperature T is related to the form of the electron velo- 
city distribution function f(v) r for the case when the plasma is immersed in 
a magnetic field by 


T 

r 


1 

k 


/- 

o 

J v 2 (v) + 


v(v ) 

2 (v) + (<D - O^) 2 


) 


(v) + (<n- o^y 


f (v ) v dv 


df(v) v dv 
d u 


where k is the Boltzmann constant, v the electron velocity, v (v) the electron 
collision frequency with neutral particles, u the electron energy, co the radiation 
frequency and the electron cyclotron frequency. 

In order to use this expression it is necessary to assume a form for f(v). 

The Maxwellian velocity distribution is of the form f(v) ~ exp(-mv 2 /2kT) so that 
as 3 . simple form of a more general d.is’tritju'tion function we will use 
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f (v) ~ exp(-bv a ) 


We then find that b 


is related to the average energy <u> by the expression 


b 


a / 2 


1 

2 


m T (5/ a ) 

<u> r ( 3 /0 


where T(5 /a ), T(3/a ) are non-integer Gamma functions. The difference 

between various forms of the distribution function is illustrated in Figure 


For nitrogen the electron- neutral collision frequency v (v) is an increasing 
function of v over almost all electron energies of interest. When T r is calcu- 
lated it is found to peak at = cu for a distribution represented by a ) 2 
remain constant for a = 2, and show a dip at cu = co for a ( 2. 


Figure 5 shows how the exponential parameter a and the average energy < u > 
changes as a function of time in the nitrogen afterglow. These results were 
obtained by matching the experimental results' against the theoretical curves. 

The time required for the parameter a to return to a value of 2, corresponding 
to a Maxwellian distribution is over 20 microseconds, corresponding to more than 
1CT electron- neutral mean free times. The electron temperature is still an 
order of magnitude higher than that of the neutral gas when a Maxwellian distri- 
bution is reached. 



Figure 4. Various calculated 
forms of the density distribu- 
tion function v^f(v) versus 
electron energy. The values of 
the average electron energy (u^> 
is related to the parameter a 
in the form chosen to represent 
f(v) ~ exp(-b v° ). a = 2 and 4 
correspond respectively to Max- 
wellian and Druyvestyn distribu- 
tions. a = 1 represents a dis- 
tribution in which more electrons 
have energies above the average 
than is the case for a Maxwelliar 
distribution. 
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Figure 5. Experimentally determined, variation in the mean energy <u> and 
exponential parameter a as function of time in the nitrogen afterglow. 


ELECTROMAGNETIC WAVE PROPAGATION THROUGH MAGNETOPLASMAS - THEORY 


An original theoretical treatment of polarization transformations in terms 
of a matrix representation Eased on the concept of orthogonal polarizations has 
been developed. Part of this was reported in the Proceedings of the IRE (Refs. 

18 and 19) and a fuller treatment was given in Technical Report 8 of the 
Rensselaer Plasma Laboratory (Ref. 20). 

Relations were developed for propagation properties of an electromagnetic 
wave along the magnetic field direction, both for a uniform and an inhomogeneous 
plasma. Techniques were described for predicting the wave polarization which 
propagates unchanged through the medium and the wave polarization which propagates 
with maximum intensity. The expressions derived are briefly summarized below: 


The Transformation Matrix 


Since only 'two possible modes corresponding to two possible polarizations 
can propagate in an anisotropic plasma, then a wave with any other polarization 
incident on such a plasma will be split into two parts with independent pro- 
pagation characteristics. Thus the polarization of the incident wave will be 
continuously transformed as it propagates through the plasma due to the fact 
that the refractive indices of the two modes are different. We shall assume 
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that a transformation matrix M exists which transforms the input polarization 
I it > into the output polarization | jt'> . That is, 

I jt ' > = M | jt> 

The assumption of the linear behavior of the medium is generally valid if 
the propagating wave does not alter the characteristics of the medium. It is 
therefore a "small signal" theory which could be expected to be valid for 
communication problems where the energy density of the wave is much less than 
the energy density of the medium. 

We shall first investigate the theoretical form which the transformation 
matrix M assumes for propagation along the magnetic field in an anisotropic 
plasma, then present a method for measuring the transformation matrix and dis- 
cuss the propagation properties which can be determined from such a measurement. 


Propagation Along the Magnetic Field 


If a wave of arbitrary polarization is propagating along the magnetic 
field in a uniform anisotropic plasma, how does the polarization of the wave 
change as the wave propagates? 

There are two ways of answering this question depending on whether we 
describe the polarization in terms of the polarization parameters or the 
elliptic parameters. We will consider the description in terms of the 
polarization parameters. 

An arbitrary polarization can be represented as 


| Jt > = COS (9/2) 


1 

tan(©/2)e^ 


We can let | jt > represent an electric field vector of unit intensity. As the 
wave propagates along the magnetic field the polarization | jt > will change and 
the resultant polarization at any point |jt'> is determined from 

I jt ' > = M | jt > 

and our objective is to find the form of M. 

Now the only modes which can propagate along the magnetic field in an 
anisotropic plasma are the right and left circular waves, corresponding to 
two of the eigenvectors. Thus, in the two-dimensional principal coordinate 
system, we can write the polarization of any wave as 
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where the components C r and C are given by 


C r = (l//2)[cos(©/2) + sin(o/2) sin 0 - i sin(©/2) cos 0} 

and 

C t = (l/'/~2){cos(©/2) - sin(e/2) sin 0 + i sin(©/2) cos 0} 


propagates with a wave number 


k = — n 
+ c + 


! /■* 

while C t propagates with a wave number 
k =- n =- /T 

P _ (7 V 


Therefore, at any point along the magnetic field the right and left circ ular 
components of the wave will be given by 


and 


= C 


ik_x_ 

e ~ 3 


We can write these equations in the form 




* 


- 1 • 

c 

r 


ik x, 
e + 3 

0 

i 

C* 


ik x _ 
e - 3 

ik x„ 
e + 3 

- «■ 





or 



I jc' > = M 1 ~ > 


( 1 ) 

( 2 ) 
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and 


Now the polarizations | jt )> and are related by 

I * > = U |it) 


I*' > = u | *’> 

where U is the two-dimensional unitary transformation given by 


U = 


/T 


i 

i 


-2 

2 


The inverse of U is given by 


U 


-1 


/T 


2 -2 

If we substitute (3) and (1)-) into (2) we obtain 
U | jt' > = M U |jt> 

Premultiplying by U ^ we obtain 

I jt' > = U' 1 MU | jt> 

By comparing ( 7 ) with ( 5 ) we see that 

M = U" 1 M U 

If we now substitute M from (l) together with ( 5 ) and (6) we find that 


M = i 


ik x_ ik x ik x_ ik x_ 

(e + 3 + e " 3 ) -i(e + 3 - e “ 3 ) 


ik x ik x, 
i(e + 3 - e " 3 ) (e 


ik x_ ik x_ 
+ 3 H- e - 3 : 


This is the result which we sought. Its use in (2) will enable us to des- 
cribe the change in polarization of a propagating wave. 



To generalize the solution of the propagation problem in magnetoplasmas 
we must take into account the fact that the plasma is probably not uniform 
and not infinite in extent and that propagation may occur in some arbitrary 
direction relative to the magnetic field direction. We will show that the 
information obtained from two polarization measurements after transmission 
through a magnetoplasma is sufficient to answer the following three questions: 

(1) How will the plasma affect the propagation of a wave of any 
polarization? 

(2) What polarization will propagate through the plasma without having 
its polarization altered? 

(3) For what polarization will the wave leaving the plasma have the 
greatest intensity? 

This can be done in terms of orthogonal polarization measurements. 

The transformation matrix M characterizes the anisotropic medium. Thus 
the polarization | it y which is orthogonal to | n y will transform according to 


I = M 1 jr > 


(8) 


The output polarizations 


I n 1 ) and ln'> will not, in general, be orthogonal. 


If we post- multiply (2) by 
equations, and use the relation 
we obtain 


I jt y and (8) by | it > 

I I + 1 * 


|3t' >< Jt| + |Jt'> < * | =M 


add the resulting 
= 1 


(9) 


which is an expression for the transformation matrix M in terms of the two out- 
put polarizations |«' > and | «' > which result from propagating the two 
orthogonal polarizations < it l and < ~ | through the anisotropic medium. 

The transformation matrix M can therefore be determined for any given medium 
from measurements of the two output polarizations | *’> and | It ' y , to- 
gether with measurement of the intensity and phase shift changes of the trans- 
mitted signal. Once the transformation matrix has been determined for a given 
medium one then knows how any other polarization will transform. 


The Characteristic Polarizations 


From a measurement of the transformation matrix it is possible to determine 
the polarization of those waves which will propagate through the medium without 
having their polarization altered. Such a polarization will be called a char- 
acteristic polarization and will be designated by I C > . This characteristic 
polarization can be found by the usual method for determining principal axes. 

For this case the output polarization I C 1 ^ is simply some scalar X times 
the input polarization | C ) . That is 
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I c* > = M)C> 


X I C > 


from which 

(M - X 1) I C > =0 (10) 

which is a set of two homogeneous equations. A non- trivial solution exists 
only if 

Det (M - Xl) = 0 (11) 

Having measured the transformation matrix M the quadratic equation (ll) 
can then he solved to give two values for X . Using these values of X in (10) 
the characteristic polarizations I C^> can then he determined. 


The Polarization for Maximum Transmission 


The intensity of the input wave will be defined by 

I = < Jt | it > (12) 

Note that as defined here I is the sum of the squares of the vertical and 
horizontal components of the wave. If we use the normalized polarizations, 
then the input wave will be of unit intensity. 

The intensity of the output wave is given by 

I' = < 7T' I JO (13) 

We wish to find that input polarization which will result in the maximum out- 
put intensity. Using (2) in (13) we can write the output intensity as 

I' = < it | M + M | ic > (14) 

where M + is the Hermit^an conjugate of M. Note from (l4) that if M is 
unitary (i.e. if M+ = M~ ) then the output intensity I' is equal to the input 
intensity I given by (12). Also note that since the complex conjugate of I' 
is given by 


* 

I' 


+ •* 

^ it | M M i jt ^ 

< it | M + M | jt > =1’ 


then I' is real. 
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We can write the polarization ratio as 


Q = tan (e/2)e^ = r + iq 
We can write the input polarization as 


11 > 


/ l+r 2 +q 2 


r + iq 


( 15 ) 


Using (15) and the corresponding expression for | jt y we can write the 
output intensity from (l^) as 


I' = 


, 2 2 
1 + r +q 


j^l r-iqj M' M 


r+iq 


(16) 


From the measured transformation matrix M we can calculate I' from (l6). 
The quantity I' will he a real function of r and q. That is, 

I' = I’(r, q) 

The maximum value of 1 1 will occur when 


dl' 

c>r 


0 and 


Sr 

ai 


0 


By calculating those values of r and q which satisfy these conditions 
the polarization which will result in the maximum output intensity can then he 
determined from (15). The distinction between a maximum and a minimum can be 
ascertained by either investigating the second derivatives or by using the M 
matrix to actually calculate the intensity of the output wave. 


ELECTROMAGNETIC WAVE PROPAGATION THROUGH MAGNETOPLASMAS-EXPERIMENT 


Free Space Propagation 


Free space propagation at K-band frequencies through an argon plasma in a 
longitudinal magnetic field was used to verify the theory developed in the pre- 
ceeding section. Construction details of the magnet used were reported in an 
article in the Review of Scientific Instruments (Ref. 2l). The transmitting 
and receiving antenna were a pair of conical microwave horns and the polarization 
of the transmitted signal could be varied by means of the polarizer. This con- 
sisted of a length of circular waveguide whose cross-section could be varied from 
circular to oval and which could be rotated relative to the rectangular input 
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waveguide. The polarization emerging was measured using a trimode turnstile 
junction (Ref. 22).. 


Polarization measurements at 20 Gc/s were carried out for an argon plasma 
at a pressure of 260 microns immersed in a magnetic field of 4500 gauss. 


The output polarizations which result from horizontal and vertical input 
polarizations are shown in Figure 6. From this figure we find that for the 
initial horizontal polarization 


I 51 > 



and 


1.06 

0.31 /-27° 




Figure 6 Transformation of orthogonal 
polarizations. Case (a) initial polar- 
ization is horizontal. Case (b) initial 
polarization is vertical. 
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while for the initial vertical polarization 


and 



l*’> 


0.23 / 150 ° ' 
1.03 


It is important to note that the same relative phase relation exists 
between the E^ component of i it > and the Eg component of | jt'^> as 

between the E2 component of I IT and the Eg component of I H ' y 

This fact cannot be obtained directly from Figure 6 , but can be determined 
independently by comparing the phases of the signals representing the E. 
component of 1 n' ) and the E2 component of l 1 with the phase or 
a reference signal. 

Using these results in (9) we can readily calculate the M matrix to be 


M = 


1.06 

0.23 7150 

0.31 

/-27 0 1.03 


( 17 ) 


The characteristic polarizations can now be determined. The two eigen- 
values are 


= 0.92 - i 0.23 

and 

= 1.17 + i 0.23 

The characteristic polarization corresponding to is given by 

Eg/E 1 = 1.20 /89 0 = tan (e/2) e 1 ^ 

which corresponds to a right-handed wave whose orientation angle is 86 ° and 
whose axial ratio is . 84 . 

The result of propagating this polarization through the plasma is shown in 
Figure 7 a. The two polarizations are superimposed on each other indicating that 
the initial polarization was unaltered. 
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Figure 7(a) Transformation of the characteristic polarization. Input and 
output ellipses coincide. 


I |MI 

Is; 


mmm* 

IB 


mm' 
\mm i 






Figure 7(b) Transmission of the appropriate input polarization (inner 

ellipse) which yields minimum output signal intensity (outer 
ellipse ) . 

Figure 7(c) Transmission of the appropriate input polarization (inner 

ellipse) which yields maximum output signal intensity (outer 
ellipse ) . 
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The polarization which results in the maximum transmission can also be 
determined. Using (17) and (l6) we can write the output intensity as 


I' = 


1.22 + l.llr 2 + l.llq 2 + 0.l4r - -54q 

2 2 
1+r + .q. 


The condition ^ 1 1 / c) r = 0 yields 

q 2 - r 2 + 7- 7rq - 1. 57r + 1 = 0 
while the condition Si ’/Sq = 0 yields 

q 2 - r 2 - 0 . 52 rq - .4lq -1=0 

These equations can be solved simultaneously for r and q. The two 
solutions are 

r = 0 . 16 , q = - 0.82 


and 


r = -0.27, q = 1-19 

From (15) we find that these two solutions correspond to the polarizations 

0.77 

0.84 /-79° 

and 

0.87 

0.83 /-79' 

respectively. The polarization given by (l8) results in a minimum intensity 
while ( 19 ) results in a maximum intensity as can be seen by applying the M 
matrix (l7) to these cases. 

The results of propagating waves with the polarizations (l8) and ( 19 ) 
through the plasma are shown in Figure 7b and 7c. 


I*> 


max 


( 18 ) 


(19) 


I it > . 

' min 
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Guided Wave Propagation 


An X-band all metal waveguide cell was developed which incorporates modi- 
fied turnstile junctions, allowing measurement of the change of polarization of 
microwaves propagating through a magnetoplasma. The work was reported in an 
article in the Review of Scientific Instruments (Ref. 3). 

This metal waveguide cell has advantages over the type of discharge vessel 
used by other workers for propagation studies in longitudinal magnetic, fields. 
Glass discharge vessels enclosed in metal waveguide have several disadvantages: 
construction and assembly are more difficult than with an all-metal system, and 
uncertainties in the appropriate boundary conditions are introduced by the 
presence of the glass. Where a longitudinal magnetic field is present, it is 
also advantageous to have both discharge electrodes parallel to the axis of the 
vessel and to the direction of microwave propagation. This axial electrode 
symmetry minimizes wall effects. It has been found experimentally that with 
asymmetrical electrode construction, similar to that used by other workers, the 
plasma may be displaced in the magnetic field to the extent of striking one of 
the walls of the glass vessel. 

The signals from the four waveguide output ports (see Pig. 8) provide infor- 
mation about the change in polarization of the propagated EM wave due to the 
magnetoplasma when the discharge cell is operated in a longitudinal magnetic 
field. The response of crystals mounted on a pair of orthogonal ports at the 
output end can be used to obtain the magnitude of the horizontal and vertical 
components. The outputs of the other two ports are combined at a hybrid junction 
with the difference in electrical length between the two ports and the hybrid 
junction adjusted to a multiple of jt/2. The magnitudes of the signals from 
crystals mounted on the H and E arms of the hybrid junction give the magni- 
tude of the right and left circular components of the wave. 

If the amplitudes of the circular wave components entering and leaving th 
magnetoplasma are C + and C_ and C + ' and C_ ' respectively, then 

, , -(k - k_)L 

C /C = e C 

+ ' - + 



Figure 8. All metal X-band 
microwave discharge cell. 
The electrodes form part of 
the coaxial matching stub 
assembly of the modified 
turnstile junction at each 
end of the cell. 
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where 


k + = the propagation constant of the appropriate circular component 
= + 1 «+ 

P + = the phase shift per unit length 
CL + = the attenuation per unit length 
L = the length of the magnetoplasma in the direction of propagation. 

It is possible to evaluate the phase shift and attenuation parameters separately. 
The rotation of the plane of polarization per unit length is given by 


*= (P + - P_)/2 


When the electron collision frequency v is very much less than the microwave 
angular frequency 03, can be used to determine the average electron density 

n e directly. It is necessary to know the radius of the region containing the 
plasma relative to the radius of the guide, so that visual inspection of the 
plasma column in the magnetic field is required. 

When 03 is close to the cyclotron frequency, \p will be dependent on both 
n e and v . The collision frequency v may be calculated from the attenuation 
factor. 


Microwave Lens System 


A microwave lens system operating at K band frequencies was developed for 
plasma diagnostics. This work was reported in a MSA report (Ref. 23). Two 
types of lens were considered, the dielectric lens and the Fresnel half-period 
zone plate (Ref. 2k ) . The superiority of the Fresnel plate was demonstrated in 
better control of bandwidth and gain and in concentrating most of the microwave 
energy close to the axis through the center of the lens. 


PLASMA STABILITY AND TRANSPORT IN MAGNETIC FIELDS 


Analysis of the Modes of the Helical Instability of 
the Helium Positive Column 


Experimental investigations have been carried out on factors affecting the 
stability and containment of a plasma in a longitudinal magnetic field. The 
work has been done on positive column discharges in helium at pressures of ^ 1 
torr in which, above a certain critical value of the magnetic field B c , a helical 
instability develops. A mode analysis of the instabilities generated has been 
made and the effects of a small transverse magnetic field studied in detail. 

This work has been reported in an AIAA paper (Ref. 25 ) and in Physical Review 
Letters (Ref. 26 ). 
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Two sets of four radial probes, each set being in a plane perpendicular to 
the current flow, were used. The tn = 0, 1, and 2 modes were identified by 
appropriate phase shift measurements. 

The m = 0 oscillation vanished when the applied electric and magnetic fields 
were precisely aligned but could be generated if a transverse field component as 
small as 0.5 gauss was present in a total magnetic field of around 2 kilogauss. 

The amplitude of the m = 0 and m = 1 mode oscillations, shown in Figure 9 as a 
function of the transverse magnetic field demonstrate clearly the effect of B^ 


B r (gauss) 



on the m = 0 mode. The amplitude of the two oscillations is similar provided B^ 
is of order 10 gauss. It is also found that in the region outside the longitud- 
inal magnetic field the amplitude of the m = 1 oscillation dies away whereas that 
of the m = 0 oscillation remains almost constant. 

The observed m = 0 oscillation is interpreted not as an unstable mode of the 
helix as in the m = 1 oscillation but rather as a nonlinear effect arising from 
the interaction of the helical mode with the asymmetry introduced by the trans- 
verse magnetic field. 

The effect of Bj on the stability of the m = 1 and m = 2 modes was analysed 
mathematically and a modified form of the Kadomtsev criterion obtained which pre- 
dicted mode switching for a given longitudinal field as B T was increased. The 
experimental results shown in Figure 10 verify this analysis. The effect of 
small transverse magnetic fields on the overall stability considerations of the 
plasma has thus been shown to be quite significant. 



Figure 9- Upper: amplitudes of 

the oscillations in discharge 
current (i ) and the m = 0 
and m = 1 modes as functions 
of the magnetic field. Lower: 
coupling coefficients, 

C Q = |m=0| / |m = lj and 

C.=I / )m = l/ ,as functions 

of the transverse magnetic field. 


21 



Figure 10- Comparison between theory 
and. experiment illustrating mode 
switching from m = 1 to m = 2 as 
the angle of inclination 9 between 
the discharge tube and the magnetic 
field is increased. The points 
represent experimental data. The 
solid line is the theoretical curve 
for a pressure of 0.3 Torr, normalized 
at 9 = 0°. 


© (degrees) 


Subthreshold Excitation 


It is common practice to use a linearized theory to predict the critical 
magnetic field B and the characteristic of the fully developed helical in- 
stability, even though nonlinear effects must be involved in the anomalous state 
of the positive column which is studied experimentally. A subthreshold experi- 
ment was carried out by exciting a pair of probes in the discharge and measuring 
the received signal at another pair of probes 40 cm farther down the column as 
the magnetic field was increased to a value just below B^. 

A typical result is shown in Figure 11. A well-developed resonance in the 
received signal occurs as B nears B c - The sub- threshold resonance is within 1 tfo 
of that of the self-sustaining oscillations which the plasma exhibits at B = B c . 
This experimental check serves as a test of the linear assumption in the 
Kadomtsev theory (Ref. 27). This work has been reported in an Applied Physics 
Letter (Ref. 28). 

The sub- threshold technique was also applied to a study of the separate 
m = 1 and m = 2 modes of the Kadomtsev instability. The theoretical formulation 
derived previously indicates that when the transverse field component B T exceeds 
a certain value the m = 2 mode becomes less stable than the m = 1 mode. Experi- 
mentally obtained sub- threshold response spectra confirming the effect are shown 
in Figure 12, indicating that the m = 2 mode is a genuinely independent mode of 
instability. 
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Signal 



B : 1500 gauss 



Figure 11 . Signal received at second 
probe set as a function of the fre- 
quency of the source connected across 
the first probe set. At the critical 
field strength (2100 gauss), the 
plasma exhibited self sustaining 
oscillations at 16.3 Kc/s. 



Received 




Figure 12 . Sub- threshold response 
spectra as a function of the applied 
transverse magnetic field. The arrows 
indicate the frequency of self-sus- 
taining oscillations that appear when 
the magnetic field is increased above 
the critical value (B ). 
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Minimum B Geometry 


Preliminary experiments aimed at stabilizing the plasma by means of quadrupole 
magnetic fields (an Ioffe bar configuration) have shown a marked effect as indi- 
cated by Figure 13. However this geometry is difficult to analyse mathematically. 


w . o / Figure 13- Stabilization of the dis- 

v probe- *ycm charge by a quadrupole magnetic field. 

Top trace: AC voltage picked up by a 

l s tab ■' 40 A/cm probe Bottom trace: 

current through the quadrupole windings. 



5 ms/cm 


A simpler magnetic field geometry which should, in principle, yield significant 
results in terms of the development of the helical instability is that of an 
azimuthal B field. Work is in progress on setting up such a system using a central 
current-carrying conductor along the axis of a hollow annular discharge cell. 


Cesium Plasma Cell 


A special cesium plasma cell has been constructed which has several novel 
design features. There are three distinct regions each characterized by separate 
vacuum systems: two electron gun assemblies and a neutral plasma region. The 

plasma source consists of a cesium reservoir with 1000 hollow stainless steel 
rods which inject the neutral cesium atoms directly on to the hot tungsten col- 
lector. When run double-ended, that is, with cesium sources at both ends of the 
cell, in a magnetic field of order several kilogauss, it is expected to achieve 


densities of the order of 10 


14 


cm 


-3 


The Initial Value Problem of Plasma Transport Phenomena 


A theoretical study has been made of the initial value problem involved in 
the mathematical treatment of rapidly varying plasma phenomena. The transport 
of electrons from a given initial state in a neutral- collision- dominated plasma 
is treated for the cases of diffusion, flow in a uniform electric field, and 
flow across a magnetic field. This work has been reported in a MSA Contractor 
report (Ref. 29). 
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The adequacy of the classical theory of diffusion derived from Fick's law, 
(particle flow proportional to density gradient), has keen questioned when it 
is applied to plasma phenomena which vary rapidly with time. The macroscopic 
telegrapher's equation includes the effect of particle inertia but does not 
necessarily provide a more accurate representation of the particle flow than 
Fick's law. It does however show that flow of fluid from a region which is 
initially at a higher density to other regions involves a particle flow which 
decays exponentially in time, and that the flow has a leading edge that moves 
with the thermal velocity. There is no prediction of an infinite velocity of 
propagation as in the formulation derived from Fick's equation. 

As illustrated in Figure 1^ the telegrapher's equation describes a flux of 
particles moving at the thermal velocity which is not predicted by the classical 
equation. While this appears to be a refinement in the description it also serves 
to point up the defect of all macroscopic descriptions because the leading edge, 
which apparently is so sharp, is in fact drawn out in space in a way which depends 
upon the number of particles in each velocity class in the initial distribution. 

To recover an accurate description of the particle flux we must therefore resort 
to a microscopic formulation of the problem with the velocity distribution function 
explicit and with particular attention paid to the initial conditions. We use 
the Boltzmann equation to discuss the problem of time- dependent diffusion from 
given initial conditions in a Lorentz gas in the presence of external fields and 
spatial gradients. 



Figure 1^. Spatial diffusion of particles from an initial rectangular diffusion 
in configuration space. Comparison between the Fickian (upper) and telegrapher's 
equation (lower) descriptions is shown. 
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The Boltzmann equation for the electrons can be written in the form: 


df 

at 


+ v 


af 

a r 


+. R 


af 

a^ 


/ av » f 


do s(x)g(f ’F'-fP) 


where R is the external force per unit mass, F(r, v) the number density of 
neutral particles on the space df (r, v) which is assumed to be constant in 
time, the neutral particle velocity, g the relative speed between the 
electrons and neutral particles, s(x) the scattering cross section in terms 
of the scattering angle X and do = sin X dX d0 is the element of volume in 
scattering space. The primed functions are evaluated at velocities that are 
pre-collision with respect to v and Vjj which are then regarded as post-col- 
lision velocities. Since we ignore the recoil of the massive neutral particle, 
F is equal to F', g = v and v' depends only on X. Then the v^ integration 
can be replaced by the density N of neutral particles in physical space. 


For simplicity we restrict the treatment to the case where R and N do 
not vary in the spatial y and., z directions. Then if f satisfies the 
Boltzmann equation, so do the y and z derivatives. We take these deri- 
vatives to be zero every where initially, and it follows that they remain zero. 


The integro-differential equation may be expressed as a set of differential 
equations by using spherical coordinates (v, 9, 0) in velocity space and 
expanding f in spherical harmonics 


f(x, v, t) = 



z 

£ f |» <*• r j!„ <®> <»■ 


A hierarchy of equations is then obtained in the f^| m . ■ For the description of 
deviations from isotropy that develop in a limited time it is appropriate to 
focus attention on f QO (denoted by f D ) and the fp m , neglecting terms of higher 
order in Z . Introducing the real functions 


f f H + f l, - 1 

1 V2 


f 2 = 


f - f 
11 1, - 1 

V2 i 


and 


f = f 
3 10 


we obtain a truncated hierarchy of equations. 

In the absence of a magnetic field, it is convenient to use the Boltz- 
mann equation in Cartesian velocity coordinates, v = (vq, V 2 , V 3 ). Then f is 
f 0 + v"l (v.f) where f is (fn, f2, fg). If f is initially isotropic then we 
always have f = g Q + g 1 where g Q = f is the part of f that is even with 
respect to v^ and gq = v -1 (vq f-|_) is the part of f which is odd. 

The truncated hierarchy then becomes 
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Diffusion 


The application of this approach to free diffusion will he briefly indi- 
cated in this report. Further details are given in reference 30 where the 
case involving electric fields is discussed. 

In the case of diffusion £ = 0 and the two equations in the hierarchy may 
he combined to give 
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a s i 

bt 2 
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a Si 



iii 

a x 2 


for both g and and hence for f itself in the present approximation. It 
is interesting to note that this is a telegrapher's type equation at the micro- 
scopic level. 

For an electron fluid composed of n Q electrons initially concentrated in a 
delta function so far as the x-variations are concerned, the initial density 
and its time derivative are 

n (x,0) = n , 5 (x) and n (x,0) = 0 
O u 

where n^ is the time derivative of the density. Further, we use the initial 
condition of an isotropic distribution in velocity space 

g Q = F(x,v) and g;L = 0, 

which imply that 

f (x,v,0) = F and (-f§ ) = q = ' v i "H 

Also F(x,v) = n Q 5 (x) G(v) 

where G(v) is the initial velocity distribution. 
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It follows that f = 0 for |x| > | v | t and elsewhere f is given by 


f - .xp(--|S)G(v)n o 5 (,-vt). » 

where the Bessel functions are evaluated at 

use the modified Bessel functions which are 

functions by I (x) = i n J (ix). 

J n n v ' 



. . _i -/ 

["/(S) v-l) ^vf t 2 - x 2 J • We 

related to the ordinary Bessel 


This equation describes a fluid composed of electrons of a particular 
velocity class which shoots out a leading edge followed by a residue that be- 
comes Fickian after long times. Thus the behavior is similar to that described 
by the macroscopic telegrapher's equation. On the other hand the macroscopic 
density derived from the two approaches can be substantially different. For 
example, if the initial distribution of velocities is Maxwellian, then the 
density of particles associated with the free- streaming is such that there is 
no leading edge. This is illustrated in Figure 15 which points up the differ- 
ences between the two macroscopic approaches and the microscopic approach in a 
qualitative way. The separation of collided and uncollided particles is made 
explicit. 


It is important to notice that the description based on the macroscopic 
telegrapher's equation is least accurate at early times where it differs most 
from the classical description. In fact for a given initial condition there 
may be no time at which the macroscopic telegrapher's equation adequately 
represents the correction to the classical diffusion equation. It is not 
possible therefore to identify any version of the macroscopic theory as the 
preferred version during the early stages of diffusion where the descriptions 
differ. The flow depends strongly on the initial distribution of velocities, 
which is hidden in the macroscopic descriptions. Moreover the constraints used 
to terminate the hierarchy of macroscopic equations, which must be arbitrary to 
a certain extent, can have a significant effect upon the result, particularly 
for short times. In order to derive an appropriate description we conclude that 
a microscopic analysis must be followed. 
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Figure 15- Temporal diffusion of particles 
past a point x for a gas initially 
localized at x = 0. Comparison between 

(a) the classical or Fickian description, 

(b) the macroscopic description based on 
the telegrapher's equation and (c) the 
micro.scopic description based on the 
Boltzmann equation with a Maxwellian 
distribution of initial velocities. 
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ATTEMPT TO GENERATE PLASMA OSCILLATIONS WITH A 
RELATIVISTIC ELECTRON BEAM 


It has been known for many years that plasma electron oscillations occurring- 
in gas discharge vessels can be detected with pfobes. More recently, it has been 
shown that these oscillations can cause radiation that can be detected by a sen- 
sitive receiver placed outside the discharge vessel (Refs. 3Q>31 )• Looney and 
Brown (Ref. 32) showed that, at least in their specially constructed vessel, the 
frequency and amplitude of plasma oscillations depend strongly upon the state of 
electrode sheaths. However, any modulated beam, passing through a plasma, may 
set up plasma oscillations, providing one of the Fourier components of the beam 
is an rf current whose frequency is close to the plasma frequency. Such modulated 
beams would be expected to exist in discharges because of the beam- plasma insta- 
bility (Refs. 33,34,35). 

An experiment has been performed in which a bunched beam from the Rensselaer 
Linear Accelerator was fired through a decaying plasma. A sensitive receiver was 
placed outside the vessel to monitor the level of microwave radiation resulting 
from the beam-plasma interaction. 

This experiment possessed the following features: 

1. The use of a decaying plasma eliminated the possibility that external power 
circuits were furnishing the power to maintain oscillations. All the power for 
oscillations had to come from the beam. 

2. Since the plasma was allowed to decay, it approximated a cold plasma. 

3. By the very nature of the accelerator, the beam was pre- modulated, so that 
sheath modulation or instability modulation by the plasma was not required. Con- 
versely, the beam electrons moved with almost the speed of light, so. that any 
oscillating fields in the plasma could not impress additional modulation of the 
beam density. Furthermore, the beam was not injected through an energized elec- 
trode, as in the Looney-Brown experiment, or other experiments using a separate 
electron gun. 

4. The beam could be injected, into the plasma at any pre-selected time during 
the plasma decay in 0.1 microsec. bursts. During this time the plasma density 
decayed by about l/lOO of 1 ’jo. Also, because the high energy electrons have a 
very slight tendency to ionize the argon gas, the beam did not reverse the decay 
of the plasma density. 

5. Probes were not inserted into the beam- plasma interaction region, since it 
was feared that these might perturb the plasma or account for spurious detection 
of oscillations. 

The sixteenth harmonic of the beam's fundamental frequency is 20.800 kMc. 

The Linac beam was fired through the decaying plasma when the plasma frequency 
near the center of the vessel was close to this frequency. 'If the bunched beam 
from the accelerator is regarded as the Fourier series superposition of a dc 
current and rf currents, then the amplitude of the sixteenth harmonic current 
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is about 1/3 of the amplitude of the fundamental current. It was expected that 
the current's sixteenth harmonic would interact strongly with the plasma, while 
the higher and lower harmonics would interact weakly. Then a peak beam current 
of 0.3 amp would have the same effect as a 0.1 amp rf current. 

The discharge vessel was 28 cm long and 5 cm in diameter. Metal caps at 
the ends of the vessel acted as electrodes and also as entrance and exit windows 
for the accelerator beam. The vessel was filled with argon gas to 0.115 torr 
and a discharge was initiated by a 60 microsec, 15 amp pulse from a plasma pulser. 
A microwave bridge was used to monitor the plasma density in the vessel. The 
interferometric fringe pattern, exhibited on an oscilloscope, displayed no jitter, 
indicating that the pulser created essentially the same density each time. 

The accelerator beam was fired along the axis of the vessel. The region 
of interaction between the beam and the plasma was then a cylindrical region co- 
axial with the vessel. Its length was the length of the vessel, 28 cm, and its 
diameter was the diameter of the beam, 1.6 cm (Ref. 36). 


Due to the diffusion of plasma electrons and ions to the walls of the vessel 
during the decay, the plasma electron density in the interaction region was non- 
uniform. The plasma density at the boundary of the interaction region was esti-. 
mated to be 10$ lower than the axial density, making the plasma frequency there 
5$ lower. This means that in the interaction region the beam passed through a 
distribution of plasma frequencies. 

Ten times each second the plasma was pulsed and allowed to decay. Then, at 
a pre- determined time during the decay, the Linac beam was fired for 0.1 micro- 
sec. The time delay between pulsing the plasma and firing the Linac beam could 
be varied in 2.5 microsec steps. The longer the time delay between pulsing the 
plasma and firing the Linac beam, the smaller the plasma frequency at any point 
when the beam passes through it. It was determined, by examining the interfero- 
metric pattern, that in 2. 5 microsec the plasma frequency near the center of the 
vessel changed several tens of Me. Covering the entire decay period in such 
small steps assured that the beam and the plasma would pass close to resonance. 
This was augmented by the fact that the interaction region contains a distri- 
bution of plasma frequencies. 

A K-band microwave horn was placed outside the vessel and was connected 
to a superheterohyne receiver. The receiver bandwidth was 8 Me and its sensi- 
tivity was 2. 5 x 10” ^ watts. The accelerator fundamental frequency and the 
receiver frequency were measured with wavemeters consisting of resonant cavi- 
ties. The receiver frequency was varied in 5 Me steps in the range 20,690 Me 
to 20,890 Me. 

2 

The receiving antenna had an effective area of 75 cm . It was located 
25 cm from the vessel, midway between the ends of the vessel, where it inter- 
cepted l/50 of the radiation at right angles to the axis of the vessel. The 
receiver was expected to respond to radiation at a higher level than 2*5 x 10 
watts. Using a peak current of 0.3 amp (equivalent to an rf current of 0.1 amp) 
it was not possible to stimulate oscillations at a sufficiently high level to 
excite the receiver. It can therefore be stated that the power emitted from 
the vessel was less than 1.25 x 10"° watts, or 0.4 x 10 -10 watts per cm of 
path transversed by the beam. A detailed report has been prepared (Ref. 37 )* 
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